ISU-IAP.Th93-03, Irkutsk 



7r + and 7T° polarizabilities from 77 — > tttt data 
on the base of S-matrix approach . 

A.E.Kaloshin and V.V.Serebryakov 1 

Institute of Applied Physics, Irkutsk State University, 
664003, Irkutsk, Russia, kaloshin@physdep.irkutsk.su 
1 Institute for Mathematics, 630090, Novosibirsk, Russia. 



Abstract. 

We suggest the most model-independent and simple description of the 77 — * tttt 
process near threshold in framework of S-matrix approach. The amplitudes contain 
the pion polarizabilities and rather restricted information about tttt interaction. Ap- 
plication of these formulae for description of MARK-II M and Crystal Ball [T2J data 
gives: (a - [3) c = (6.0 ± 1.2) • 10" 42 cm 3 , (a - /3) N = (-1.4 ± 2.1) • 10~ 42 cm 3 (in units 
system e 2 = 47ra) at the experimental values of tttt scattering lengths. Both values are 
compartible with current algebra predictions. 
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1. Introduction. 



The appearing data on 77 — ► tt + tt~ MARK-II M and 77 — > 7r°7r° Crystal Ball 
H near the threshold give the opportunity to obtain the pion polarizability from 77 
experiments. Recall that for the ir + polarizability there exist two measurements in 
processes on nuclei 0, U but for tt° the 77 experiments give the first possibility to 
measure it. The polarizability's measurements allow to verify the predictions of differ- 
ent low-energy models based on chiral symmetry. Note that a straightforvard using of 
chiral model for description of 77 — ► tttt cross-section does not have big sense becouse 
of huge unitary effects. One can see (see, for instance, ||) that the one-loop chiral 
model calculations || differ significally from unitary formulae and experiment even in 
region of 400-500 Mev. 

By this reason we use our S-matrix approach @ to this reaction. Note that in 
other unitary models || || [TOj polarizabilities are not controlled and there exist some 
problems with low-energy theorem even on classical level. 

The polatizabilities in the 77 — ► tttt amplitude play the same role as the electro- 
magnetic radius in a formfactor: they are the structure constants in threshold decom- 
position. So the problems at extraction of a, (3 from a cross-section are in principal the 
same as in < r 2 > extraction from pion formfactor at q 2 > 4m 2 . The main problem 
is that the polarizability's effects can not be separated from effects of tttt -interaction 
in final state and one is forced to utilize this information in the 77 — > tttt amplitude. 
Note that the pion polarizabilities can be obtained from any analysis of 77 — > tttt data 
in wide energy region by extrapolation to point s = 0,t = m 2 . Such an analysis is in 
fact rather complicated task (in our opinion it is not made correctly up to now) which 
needs first of all the better understanding of / = J — hadron dynamics. So we think 
it will be reasonable to develop the more model-independent and simple description of 
near-threshold region 77 — > tttt which includes only some well-established properties of 
7T7r-interaction. We hope that future 77 — > tttt experiments can utilize these formulae 
in physical analysis of data. 

In [JO]] we have been performed the simplest analysis of 77 — > tt°tt° [Q] data to 
extract (a — (3Y° . Here we analyse the obtained formulae and use them for description 
of data 10], 0] on both reactions. 
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2. Simplified unitary formulae for 77 — > txtx near the threshold. 

Let us introduce the helicity amplitudes T\ in c.m.s. in usial manner. The cross- 
section 77 — > 7r + 7r~ for unpolarized photons is defined as 

da p(s) ,. m ,0 , m ,2-1 

5n = 2^F {|r++l +|T+ - |} (1) 

where p(s) = (1 — 4/i 2 /s) 1//2 , ji = m^. For 77 — > 7r°7r° one should add the factor 1/2! 
in (1) to take into account the identity of 7r°'s. Isotopical amplitudes 77 — > 7r7r are: [] 

T C = T( 77 - vr+vr-) = ^- T'=° + ^ T'= 2 

T N = T( 77 - ttW) = ^ T<=° - ^| T'= 2 (2) 

It is well known that helicity amplitudes contain the kinematical singularities and 
zeroes. For this reaction we can pass over to reduced amplitudes M\ which do not 



contain them [13]. 

T ++ = sM ++ , T + _ = (tu - //)M + _ (3) 

Besides the using of amplitudes M\ is very useful in another aspect: it allows to control 
easy the low-energy theorem including the structural corrections - see below. So we 
shall use the reduced amplitudes Ma in the following. 

The starting "final state interaction" formulae for M ++ are contained in [/J. They 
satisfy the one-channel unitary condition for lowest partial wave and are normalized 
on the pion polarizability. Note that the polarizability is the attribute of just the 
lowest partial wave. The formulae of have two arbitrary constants C° and C 2 , 
if the contributions of nearest cross-exchanges are known. But one can refuse from 
this additional suggestion by introducing of some effective constants which absorb any 
high-energy contributions. Only the lightest p, u exchanges should be considered as 
" alive" . As a result we have the following expressions for reduced helicity amplitudes: 

M? + = ^(C° + HUs))-^n%s)(C 2 + HUs))+W N (s,t) (4) 



1_ We indicate these relations here to avoid the mistake as in some recent papers [10, 12] 
on this matter, where the T N in (2) has the opposite sign. The reason of this mistake may 
be found in any old textbook: you should denote \ tt + >= — | / = 1, 13 = 1 > to use the 
standard Clebsh-Gordan coefficients. 
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Here Q 1 is the Omnes function, corresponding to phase shift 5q. 

7T J S\S' — S) 

W°\s,t) = W n (s,t) + W p (s,t) + a c (5) 
W N (s,t) = W(s,t) + W u (s,t) + a N 

W" = 2g2/i2 W = -Z"( - - ) Z R = q 2 

(t — fi z )(u — fi 2 ) m 2 — tm 2 — u ' 

Functions Hy(s) in (4) are the rescattering contributions: 

a° ,a N in (5) are some effective constants from cross-contributions and possibly s- 
channel high energy contributions. If the decays widths R — > tc-/ are known, they may 
be estimated (see below) but we prefer to consider them as free parameters. 

The J=0 projection of (4) must satisfy the one-channel unitary condition for iso- 
topical amplitudes 1=0,2 . 

ImF\s) = P {s)F\s)f:i^ p(s) = (1 - A^/s) 1 ' 2 (7) 

It is easy to see that the unitary condition is reduced to 

V\s) = W^ J = (s) (8) 

where W 1 are the isotopical combinations of W c , W N (5) which are projected on J=0 
state. 

The amlitudes (4) contain four free parameters: C°,C 2 ,a° ,a N . But becouse of 
smallness of 5 2 phase shift the Omnes function Q 2 is very close to 1 and C 2 Q 2 does 
not differ from constant . So we can put C 2 = in (4). As a result we have 3 free 
parameters for two reactions 77 — > 7r°7r°. The pion polarizabilities are related 

with them as 



2 V v 7 ^ m 2 p -fi 
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^(a-(3) N = ^L + a N - Z?-^— 2 - Z«^— 2 (9) 
2 V y/3 m 2 - fi 2 ml- fi 2 

In the folllowing we use (a — (3) c , (a — (3) N , C° set of parameters. The reason of 
appearing the additional parameter C° is clear: the point of polarizability's definition 
is out of physical region of 77 — > irir and final state interaction leads to Vt°(s) 7^ f2°(0). 
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Let us note that amplitudes (4) may be rewritten in another equavalent form. If 
5q — ► 7r at s — ► oo , there exists the equality 

a^ = l + ^ + U7(^ Im cm- (10) 

and similar for 1=2 without (3s term. It leads to disappearing of constants oF , a N in 
(4), to appearing of C 2 and replacement of C° by polinom. That is an equavalent form, 
but we prefer to use the form (4). 

3. About the nix— interaction. 

Recall [0] that the formulae (4) satisfy the one-channel unitary condition for lowest 
partial wave J=0 and 5® denotes only smooth part of experimental phase shift without 
the sharp / (975) effect. For justification of such approach one should consider the 
multichannel problem (see, for instance, fI3| ) but it can be easy understood at quali- 
tative level. Simplest (and widely used) suggestion on the structure of tttt amplitude 
is (see review fl"5|j ): 

8° = 6 Back9 + 5 f °, rio = V fo (11) 

In any case the tttt — > tttt J=I=0 amplitude contains the smooth contribution with 
typical scale ~ 500 Mev and the sharp fo one with scale ~ 50 Mev. Solving the 
multichannel problem [14]] we see that these two contributions are separated from each 



other in the answer. If the coupling /o77 is very small ( that is an experimental fact, 
the decay width f — > 77 is the value ~ 0.1 Kev ) then the presence of f leads only 
to the local interference effects near 1 Gev WM and its influence is negligible in a big 



energy scale . So for region i/i < 800 Mev the effective one-channel problem is arised 
and result we have the amplitudes (4). 

So to calculate the 77 — > tttt we should use the main properties of smooth part 5® 
phase shift. Firstly that is the threshold region which is described in term of scattering 
lenghs and the chiral zeroes. 

[s - sj) (Afi 2 + gi) 
(s + Sl ) (4/i 2 - 4) 

where s% ~ m 2 . This formula was used for describing the titt data and the analysis 



I6fl gives: 

s° = -0.2/i 2 , s 2 = 2.4/i 2 , 

a° fj, = 0.24 ± 0.04, = -0.04 ± 0.04 (13) 
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Secondly, we slould utilise the large scale behaviour of phase shift. In experiment 
smooth part of 5® crosses the ti/2 at yfs ~ 900 Mev and looks as linear on y/s function 
in very wide energy region. 

In this situation we describe threshold region as (12), then the phase shift is de- 
scribed as linear on y/s function, and it reachs 7r at \fs = 1.5 Gev — see Fig. 1,2. 
Another variant of action may consist in using for amplitude im — > mr the N/D repre- 
santation type of |T7| which includes all above-mentioned properties. We can also 
modify sligtly this model for better describing of experiment. But we have been con- 
vinced that the phase shift behaviour at say 1.3 Gev practically does not influence on 
the extracted values of polarizabilities. So we use below the simplest linear approx- 
imation! with threshold behaviour (12). As for the phase shift S 2 , we need only the 
threshold properties (12) becouse it influences only on the rescattering integral which 
is saturated by the threshold region. 

4. The low— energy theorem and polarizabilities. 

It is well known (see, for instance, review [|18[]) that Compton scattering amplitude 
on hadron contains the classical charge normalization |TU| and next structural terms 



of threshold decomposition. In case of spin target there appear in uj 2 order two 
structural constants a, (3, named the electriical and magnetic polarizabilities. One can 
relate them with the helicity 77 — > amplitudes 0, which have the following form 
in the vicinity of point s = 0, t = fi 2 : 

T++ = s[Ml + + ^(a-(3r] 

T + - = (tu - /i 4 ) [ Ml_ + ^LlM. ] (14) 
Here M 71 " is the contribution of 7r-exchange (QED) for the charged pions and zero for 



neutral. Note that here the reduced helicity amplitudes m |13| are arised , which are 
useful to control the low-energy theorem. 

Up to now there exist two measurements of tt + polarizability in experiments on a 
target. The experiment nZ — > n^Z [|J gives the value (a— (3) c = (17.0±3.6)-10~ 42 cm 3 , 
the result of another experiment 7p — > / -fir + n || is (50 ± 30) ■ 10~ 42 cm 3 . 

As for different model prediction for a,/3, let us begin from current algebra (see 
review [T21]] on this matter). The using of PCAC leads to (a—/?) 71 " = 0, the polarizability 
of 7r + is reduced to formfactor tcVA, which is measured in the decay tt — > ^ev. The 
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using of last experimental data for this decay |2T| leads to: 

(a - (3f = (8 -T- 10) ■ 10" 42 cm 3 (15) 

The majority of model predictions, based on different forms of the chiral model, also 
are near this value or slightly higher. 

As for neutral pion polarizability, all models predict smaller ( or much smaller) 
value as compared with charged pion. But the value ( and sign) depends strongly on 
the model. For instance in the quark- virton model |22[ (a — f3) N = 2.8 • 10~ 42 cm 3 , 



in the four-quark model of superconductor type|23[| (a — (3) = —3.2 • 10 42 cm 3 



5. Our preliminary expectations. 

In framework of the approach |7[] one can obtain some estimates of the parameters 
in (4). The starting point of |7| is the dispersion relation for the reduced helicity 
amplitude M ++ at fixed t ~ /i 2 , which is written without subtractions according to the 
standard Regge behaviour. Putting s=0 in it and saturating the cross-integral by the 
nearest resonance exchanges we shall obtain 

2 V 3 

|( Q - pf = -Z» -Z" + Z h - + Z h - + ^ - ^-C 2 (16) 

Here C°, C 2 are the corresponding s-channel inegrals at s = 0. If to put C 2 = and 
to subtract one equation from another there arises the relation between polarizabilities 
of charged and neutral pions. But in fact in 1=2 channel there exists one contribution, 
which can not be neglected — that is the rescattering. But this contribution can be 
easy calculated becouse it is defined only by the threshold nn properties. So taking 
the experimental values (13) we shall obtain C 2 = —0.0328 Gev~ 2 and the relation 
between polarizabilities has the form: 

(a - (3) c - (a - p) N = 9.1 • 10~ 42 cm 3 (12.6 at C 2 = 0) (17) 



if to take the existing values for R — > tt-/. Q Note that the scale and sign of right 
side (17) is defined in fact by the decay width u — > which is measured with best 
2 The relation (17) demonstrates one of the possible contradictions with using the wrong 
isotopics [0, It leads roughly to the changing of sign in (17) and contradicts evidently 
to present experimental and theoretical knowledge about a, (3. 
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accuracy. As a result we can estimate (a — (3) c , C°: 

(a - p) c = 9 ■ 1(T 42 cm 3 , C° = 0.10Q Gev~ 2 , 

at {a - (3) N « {a - /3) c (18) 

We can expect such scale values for (a — j3) c , C° in case of absence of different exotic 
phenomena ( big odderon contribuion in 77 — > 7r°7r° ||25|| , some unknown dynamics in 
1=2 channel,...). Let us note also that the smallnes of tt° polarizability as it seen from 
(16) is the result of cancellation of rather big different contributions near the point 
s = 0,t = fi 2 . So it does not means the smallness of 77 — > 7r°7r° amplitude in wide 
kinematical region. 

6. The analysis of 77 — > tttt data. 

The numerical analysis of the obtained amplitudes shows that the most favourable 
region for extraction of a, [3 is \/s < 500 Mev, where practically closed description is 
arised, i.e. here the function Q°(s) is defined mainly by the threshold properties of 
phase shift. Besides the contributions of p, u exchanges is not essential here. So let us 
start from consideration of this region. 

6.1 Near— threshold region: y/s < 500 Mev . 

In this region there exist the MARK-II data flj Q ( 5 points between 350 and 400 

Mev ) and Crystal Ball [Q]. The quality of data does not allow to perform the total 

3-parameter analysis, so let us put (a — f3) N = 0. As we shall see later in all variants of 

analysis the neutral pion polarizability is compatible with zero. Besides it was found 

that the cross-section 77 — > 7r°7r° depends essentially on all three parameters whereas 

77 — > 7T + 7r~ below 400 Mev only weakly depends on the (a — [3) N . Unlike of |IT] we 

use exact value of m^o in phase volume which is essential near threshold. 

3 In the region ~ 400 Mev in the process 77 — > tt + tt~ the helicity amplitude M_| is 

not negligible. But here it practically coinside with QED contribution ^ , so we take into 

account only this contribution. At higher energy it is nesessary to construct the model for 

helicity 2 amplitude too. 

4 In amplitudes all masses are equal to . In principal the effects of mass difference could 

play some role in amplitudes too becouse of existence of QED "amplifier" in rescattering. We 

checked it on the base of K-matrix approach and found that amplitudes are sensitive to the 

exact value of 7r + mass but not to m n o . This fact is defined by the unitary constrains and by 

presence of chiral zeroes in tttt amplitudes. 
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1) . Let us try to describe MARK-II [|lj data separately. We shall obtain very broad 
distributions for x 2 , as for C° here is practically a plateau. It means that cross-section 
is sensitive only to one effective parameter — some linear combination of (a — (3) c and 
C°. This fact may be seen from Fig. 3 where the lines (almost linear) of equal \ 2 are 
depicted in plane of parameters. Any point on this line gives the same cross-section 
in region of MARK-II threshold data but at higher energies they are very different. 

2) . Describe only Crystal Ball data in this region. Again very broad distributions 
are arised, the minimum of \ 2 corresponds to (a — j3) c ~ 5 ■ 10~ 42 cm 3 . One can see 
however that the data Jl|, [2| do not contradict to each other. 

3) . Describe the data |j] and together. It leads to: 

(a - (3) c = 5.4 • 10" 42 cm 3 , \ 2 = 15.4 at NDF = 8 (19) 
The C° parameter is defind very badly - situation is similar to Fif.3. 

6.2 Fit in extended region. 

Let us extend the energy region of Crystal Ball data used in analysis up to 850 
Mev. At higher energies helicity 2 amplitude becomes essential |F| but variation of this 
bound does not change the results essentially. 

1). Fit only CB data by three parameters. 

As a result we have been found the existence of two minimums with close values of 
2 . 



X 



I. (a - (3f = 7.0 • 10- 42 cm 3 , (a - (3) N = -0.15 ■ 10- 42 cm 3 , 
C° = 0.156 Gev~ 2 , \ 2 = 9-1 at NDF = 9 



//. (a - (3f = 8.5 • 10- 42 cm 3 , (a - (3) N = 0.4 • 10" 42 cm 3 



C° = 0.026 Gev~% x = 9.4 at NDF = 9 (20) 

The second minimum leads to sharp increasing of both cross-sections above the 
region of fit which contradict to CB data here. So we shall not consider this variant 
in the following. It is interesting that both mimimums do not contradict to theshold 
MARK-II data — see Fig.3. 

2). Fit the Crystal Ball ( < 850 Mev) and MARK-II (^i < 400 Mev) data 
together. 
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Two-parameter fit, (a — (3) N = 0, gives: 

(a-pf = (6.1 ±1.2) • l(r 42 cm 3 , C° = 0.160 ± 0.012 Gev~ 2 , 

X 2 = 18.4 at NDF = 15 (21) 

Three-parameter fit: 

(a-(3) c = (6.0 ± 1.2) • 10- 42 cm 3 , (a- (3) N = (-1.4 ± 2.1) • 10" 42 cm 3 , 
C° = 0.150 ±0.020 Gev~ 2 , 

X 2 = 17.9 at NDF = 14 (22) 
The corresponding cross-sections are depicted at Fig. 4, 5, 6. 

7. Conclusion. 

So we tried in the first time to obtain an information about pion polarizabilities from 
the existing data on 77 — > ttti. Here we delibertely used the most simple description 
of amplitudes on the base of S-matrix approach. Such an approach may be called 
as "structureless" since it utilizes only some properties of the smooth tttt amplitude 
without consideration of its physical structure. We suppose that for region of < 
800 Mev such description is fully adequate and allows to extract pion polarizabilities 
from experiment. We think that more detailed investigation of this question with 
consideration of multichannel problem ( and with correction of some defects in existing 
models) will be useful. Firstly, it will allow to define more precisely the region of 
application of simplified formulae and to obtain more correct values of polarizabilities. 
Secondly, it is interesting to analyse the spectrum and properties of scalar mesons 
from joint consideration of hadron and two-photon experiments. Let us note also very 
unexpected results (from point of view of threshold characterises first of all) of CELLO 
group (12| on partial-wave analysis of 77 — > tt + tt~ above 0.8 Gev. 



As a result of our analysis of Crystal Ball and MARK-II data we can note the 
following facts: 

• The data of CB || and MARK-II 0] are described well by formulae (4) at the 
standard form of tttt interaction and do not contradict to each other. 

• The data of CB and MARK-II separately indicate on rather small value of 7r + 
polarizability which is compartible with current algebra but contradicts to exper- 
imental value (17.0 ± 3.6) ■ 10~ 42 cm 3 from ttZ — > njZ 0. The value of ~ 12.0 is 
excluded by the 77 — > nn data. 
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• At any variant of analysis the neutral pion polarizability is compartible with zero 
in accordance to current algebra. 

• At changing of scattering length = 0.24 to any side the disagreement of Crystal 
Ball and MARK-II data is arised. 

• The obtained values of a — (3 at least qualitatevely correspond to our estimates 
(18), based on the saturation by nearest resonances. So we do not see some new 
effects in such analysis. 

When this work was finished we received the paper ||, where an ettempt was made 
to extract 7r° polarizability from data on the base of multichannel model [pTOfl . Note 



that Pennington obtains (a — j3) > with much less errors than ours. But one should 
note few defects in model [f[TJ including the wrong isotopics. 
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Figure captions. 

Fif.l The phase shifts described it text at the threshold parameters (13). 
Fig. 2 The Omnes functions corresponding to Fig.l. 

Fig. 3 Lines of equal x 2 & t hting of 5 threshold points of MARK-II by two pa- 
rameters, (a — pi) N = . Different scattering lengths a[j are indicated. The error 
(one standard deviation ) coinsides approximately with distance between lines. 
Along these lines x 2 is the same: x 2 — 8-4. 

Fig. 4 Comparison of MARK-II data with best fit curve. All curves with equal 
X 2 = 8.4 coinside in this region but differ at higher energies y/s = 700 — 800 Mev. 
The helicity 2 QED contribution is subtracted from MARK-II data. 

Fig. 5 The cross-section 77 — > -k + it~ of helicity at higher energies . Two 
variants correspond to two sets of parameters (20). 

Fig. 6 The cross-section 77 — > 7r°7r°. Comparison of CB data with best fit curve 
at values (20). 
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